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We introduce an approach to the design of three-dimensional transformation optical (TO) media 
based on a generalized quasi- conf or mal mapping approach. The generalized quasi- conformal TO 
(QCTO) approach enables the design of media that can, in principle, be broadband and low-loss, 
while controlling the propagation of waves with arbitrary angles of incidence and polarization. We 
illustrate the method in the design of a three-dimensional "carpet" ground plane cloak and of a 
flattened Luneburg lens. Ray-trace studies provide a confirmation of the performance of the QCTO 
media, while also revealing the limited performance of index-only versions of these devices. 
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Transformation optics (TO) is a unique tool for the de- 
sign of complex electromagnetic media[l . TO makes use 
of the form- invar iance of Maxwell's equations to mimic 
spatial transformations using distributions of inhomoge- 
neous and anisotropic material constitutive parameters. 
TO has inspired many exotic devices, one of the most 
compelling of which being the electromagnetic "invisibil- 
ity" cloak. However, the use of TO as a design method- 
ology typically comes at high cost; the media derived 
from coordinate transformations generally involve spa- 
tial gradients in all nine components of the epsilon and 
mu tensors. Though it is possible to find a basis that di- 
agonalizes these tensors, the diagonal basis will generally 
be a function of position for all but the most simplistic 
and symmetric designs. Moreover, the required response 
is generally outside of the range of natural materials. 

Electromagnetic metamaterials (MMs) are used to ac- 
cess the extreme material parameters required by TO 
media. MMs, for example, were used to demonstrate a 
negative index of refraction [2? ] and electromagnetic 
cloaking [3^. However, the performance of these initial 
MM constructs was limited by a combination of narrow 
bandwidth and relatively large absorption. The typical 
limitation for MM designs has been the requirement of 
constitutive parameters that have a large range of val- 
ues for both permittivity and permeability. The imple- 
mentation of artificial magnetism, in particular, requires 
resonant inclusions that are inherently lossy and disper- 
sive, leading to absorption and reduced bandwidth. More 
recently, the development of coordinate transformation 
methods in optical design approaches has significantly 
advanced MM complexity: independent magnetic and 
electric responses are required in all directions for general 
TO designs, yet most MM elements provide a controlled 
response in one or two directions. Were one to attempt 
to control all of the tensor elements of a MM simulta- 
neously, multiple MM elements would need to be either 
co-located or closely positioned, introducing very com- 
plicated magneto-electric coupling difficult to control or 
even quantify using current MM retrieval techniques [4]. 

Fortunately, the enormous degree of flexibility avail- 
able in coordinate transformations can alleviate many of 



the complexities and limitations of TO media. For a 
given TO device that performs some function, there is 
typically little concern over the physical distribution of 
fields within the device itself; rather, it is the distribu- 
tion of fields external to the TO medium that is of conse- 
quence. Therefore, the freedom exists to choose any of an 
infinite number of possible coordinate transformations, 
subject to the boundary conditions, to achieve the same 
functionality while simultaneously optimizing certain de- 
sign criteria. In 2008, Li and Pendry introduced the con- 
cept of the quasi- conformal map (QCM) to T0[5 . It may 
be shown that the QCM minimizes anisotropy across a 
two-dimensional transformation. Applying QCM opti- 
mization, Li and Pendry conjectured that certain classes 
of TO media could be realized using dielectric-only ma- 
terials. The approach was applied to the design of a 
medium that would mask a deformation introduced into 
an otherwise flat, perfectly-conducting (PEC) plane. The 
QCM technique has subsequently been applied in several 
experiments, and has also been used to demonstrate that 
substantial improvements can be achieved in the redesign 
of more conventional imaging devices [6 . 

While the advantages of dielectric-only QCM-derived 
media are extremely important from an implementation 
point-of-view, such media can only control waves prop- 
agating in two-dimensions. In this paper we show that 
the QCM concept can be readily generalized for the de- 
sign of media that can control wave propagation in three- 
dimensions (3D); however, though greatly simplified, 3D 
QCM media are inherently anisotropic, due to the re- 
quirement of magnetic response in at least one direction. 
We illustrate the performance of generalized 3D QCM 
media with ray-traces of a carpet cloak and of a flat- 
tened Luneburg lens, comparing with 2D QCM media. 
The ray traces enable us to quantify the effectiveness of 
the methods using the standard metrics of optical lens 
design. Finally, we introduce a method to mitigate the 
effect of the approximations that are made with QCM 
and QCM-like mappings. 

Consider an arbitrary mapping in of the form x^ = 
X* (x*). The unprimed coordinates represent the virtual 
domain (subscript "v"), and the primed coordinates rep- 




FIG. 1. (Color Online) 3D carpet cloaks created using 2D 
quasi-conformal transformations. In (A) the mapping is per- 
formed in X and y, and z is invariant. In (B), the mapping is 
performed in p and z, while cj) is invariant. 



resent the physical domain (subscript "p"). The virtual 
domain is the domain as perceived by electric and mag- 
netic fields. For cloaks, this domain is typically vacuum. 
The physical domain is the coordinate system that ex- 
ists in physical space. The TO-derived material tensors 
(e = Jl) for this transformation are 



Det(A)- 



^A^ K] e*- 



(1) 



where A* = dx^ /dx^ is the Jacobian of the transforma- 
tion. Suppose that we wished to find a transformation 
that required no magnetic coupling. Since e* ^ = /i* -^ , 
the magnetic response can only be eliminated by the iden- 
tity transformation. However, if a transformation could 



be found such that e* ^ 



e* ^ di'j' and e^ 



1, then 



for waves polarized such that the electric field is con- 
strained to lie in the z direction, only e^ is relevant and 
the transformation could be realized with a dielectric- 
only medium. A conformal transformation accomplishes 
exactly this goal. However, a coordinate mapping be- 
tween a physical and a virtual domain can only be con- 
formal between two domains that share the same confor- 
mal module, or M = rriy/mp = 1. If M ^ 1, then the 
transformation is quasiconformal, and the in-plane ten- 
sor components of the material parameters will no longer 
be equal. Additionally, orthogonality is maintained at 
the boundary of the mapping via Neumann ( "slipping" ) 
boundary conditions (BCs). These BCs will generally in- 
troduce discontinuities at the boundaries that manifest 
as refractive aberrations and reflections. 

It is important to note that QC mappings are two 
dimensional (e. g. x and y) with the third dimension 
(z) invariant under the transformation. Therefore, the 
QCTO prescription is valid as long as the geometry is 
also invariant along z. An example of a structure for 
which a QC transform solution can be found is the car- 
pet cloak, shown in figureja). For a wave polarized along 
the z-direction and propagating in the xy plane, a dielec- 



tric only implementation can be applied; however, the 
full TO solution requires a magnetic response equal to 
the electric response in every direction. Therefore, we 
require that jj^z = e^ for the full TO medium. 

The TO approach can be fruitfully applied to systems 
such as waveguides [7 and transmission lines. However, 
many systems of great interest are invariant azimuthally; 
for such systems, it is then natural to seek an ana- 
log QCM mapping in cylindrical coordinates. Assuming 
M ?^ 1, we calculate the QCM Xp{xy,yy) and yp{xy^yy), 
and then make the substitutions Xy^p -^ py^p and yy^p -^ 
Zv,p' We then determine the material transformation 
in cylindrical coordinates [1} [HI [9]. Figure IFb) shows an 
example of such an azimuthally invariant geometry and 
map. Remembering that ^ = e, we obtain the following 
material parameters: 



60 = a//3, and e^ ~ e^ ~ /3 



(2) 



where a = |A|~^ and /3 = py/ pp. The factor /3 arises from 
the fact that the differential volume element in cylindrical 
coordinates is a function of p. The transformed material 
parameters must compensate for the dilation of space be- 
tween the virtual and physical coordinates; however, the 
transformed medium remains approximately uniaxial. A 
further simplification can be made if we are mostly inter- 
ested in a structure that controls ray propagation in the 
eikonal limit, where only the anisotropic refractive index 
is of consequence [10]. In the eikonal limit, we can thus 
make the substitutions 



0' 



2 f' 



Of, p^ 



M0 



1, and p'^ = a/ P'^ (3) 



to achieve an equivalent medium. Since most 

magnetically-coupled MMs provide a magnetic response 
in only one direction, the substitutions of ([3| ease im- 
plementation by requiring only one component of p to 
deviate from unity. The eikonal approximation may or 
may not be appropriate depending on the intended oper- 
ation. Note that neither ^ nor ([3| are isotropic material 
specifications. Therefore, a purely dielectric implementa- 
tion of this mapping cannot control waves with arbitrary 
incidence and polarization, even in the eikonal limit. 

We can understand the need for anisotropy by consid- 
ering the local dispersion relation for the system, which 
for our transformation is given by: 



kl/P' 



kl/a ■ 



kl/a -1 = 0. 



(4) 



where the wavenumber k has been normalized to the free- 
space value ko. This form of the dispersion relation is 
derivable directly from Maxwell's equations in a homo- 
geneous, index-matched (e = p) medium with cylindrical 
anisotropy. This dispersion relation can be achieved by 
an isotropic dielectric e = a only if /c^ = 0. The cylin- 
drical symmetry of the system ensures that /c^ = along 
the ray's entire trajectory, and the ray would therefore 




FIG. 2. (Color Online) Ray-tracing comparisons of (A) The 
bare PEC scatterer, (B) the index-only revolved carpet cloak 
and (C) the properly transformed (anisotropic) carpet cloak. 
The angle of incidence is 30°. The plots on the left consist 
of incident rays in the plane containing the optical axis. The 
plots on the right consist of rays off the optical axis. The 
incident rays are blue and the reflected rays are red. 



be properly transformed. This means that only rays in 
planes of constant (j) are properly directed (meridional 
rays) for a dielectric only medium. We note that the 
dispersion relation is actually fourth-order in /c, but we 
have suppressed a power of two because the two modes 
are degenerate for impedance-matched media[ll . Our 
substitution in ([3| preserves this degeneracy, and there- 
fore our prescription remains singly refracting. 

We can demonstrate the limitations of dielectric-only 
implementations, as well as the validity of 3D QCM me- 
dia via ray tracing ^llj. We begin by studying the 3D 
carpet-cloak [5] [T2HT8] , designed to effectively flatten a 
perturbation of the form z = O.lx cos(7rp)^ in the region 
p < 0.5. Ray-tracing analyses of both the dielectric-only 
and anisotropic implementations are shown in figure [] 
A ray trace from the uncloaked perturbation is shown 
for comparison. In both cloaks, Rays incident such that 
k(f) = are properly reflected at the specular angle. How- 
ever, when /c0 7^ The dielectric-only implementation 
does not properly redirect the rays. Rather, the rays 
appear to be slightly focused by the index distribution. 
Rays that do not lie on the optical axis experience a 
gradient transverse to their initial trajectories that grad- 




FIG. 3. (Color Onhne) Material parameters for the 3D QCM 
flattened Luneburg Lens. Only the tensor components that 
differ from unity are shown. 
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FIG. 4. (Color Online) (A) Ray-traces of collimated rays in- 
cident on a flattened Luneburg lens. (B) is a close-up of the 
index-only revolved Luneburg lens and (C) the properly trans- 
formed (anisotropic) lens. The angle of incidence is 30°. The 
rays are red outside of the mapped region and blue inside. 



ually bends the rays around the optical axis. Therefore, 
rays that were initially parallel come to lie on intersect- 
ing trajectories after the exit the device. Nevertheless, 
the index-only cloak clearly reduces scattering away from 
the specular direction, possibly masking the distortions 
caused by the cloak in an experimental setting. 

The anisotropic transformation described by ([2| , prop- 
erly cloaks the bump. The reflected rays exit the cloak as 
if they had reflected specularly from a flat ground plane. 
The rays do refract slightly at the boundaries of the map- 
ping as the material parameters differ slightly from unity. 
Refraction at the mapping boundary is included in our 
ray-tracing analysis, but its effect appears to be negligi- 
ble. 

The QCM technique holds great potential for its ability 
to modify and improve conventional optical devices. The 
Luneburg lens, for example, is a rotationally-symmetric, 
gradient-index medium that perfectly images two con- 
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FIG. 5. (Color Online) Ray- intercepts of the focal planes for 
three flattened Luneburg designs. 50 collimated rays illumi- 
nating 90% of the aperture are traced through the lens. The 
angle of incidence is 30° . A scale bar is included to show the 
spot size in terms of the unflattened lens diameter, D. 



centric spherical surfaces on to one another [W. If one 
of these surfaces is taken to infinity, then parallel rays 
are focused to the finite surface. The permittivity distri- 
bution for the Luneburg is given by [19^ ei = 2 — (r/a)'^^ 
where a is the radius of the lens. 

Despite its powerful imaging capabilities, the Luneb- 
urg is seldom used since the spherical focal surface is 
incompatible with conventional planar detector arrays. 
It was proposed that TO could be used to flatten the fo- 
cal surface and retain the same imaging performance [F. 
The flattened Luneburg was subsequently realized by us- 
ing the QCM method in 2D [6]. A 3D version of the lens 
can be obtained, following the procedure given above, by 
a medium whose non-unity parameters are given by: 



Cp = 4 = e^/^^ ^0 = ^i^^ and /i'^ = a//3^ 



(5) 



where it is understood that ei = ei{py^Zy). We consider 
a lens that has been flattened to provide a field-of-view 
(FOV) of 90°. The reduced material parameters for this 
QCM are shown in figure [] We also examined an index- 
only implementation by making the assignment e^ ^ e. 

Ray-traces depicting the focal properties of both lenses 
are shown in figure [j The collimated rays are incident at 
an angle of 30°. The anisotropic lens in figure []^b) shows 
good focusing at 30°, whereas the index-only lens in Jc) 
shows severe off-normal astigmatism. This is clearly seen 
in a plot of the ray intercepts at the flattened focal plane 
(figure [] Rays parallel (but not on) the plane contain- 
ing the optical axis are over-steered such that they cross 
the plane containing the optical axis before they inter- 
sect the bottom of the lens. Even though the anisotropic 
lens offers better performance than its index-only coun- 
terpart, there are still aberrations present due to the uni- 
axial approximation in|3] This approximation effectively 
stretches the virtual domain in one direction [1 7 , such 
that the Luneburg lens is no longer spherical in the vir- 
tual domain. 

For a FOV of 180°, the QCM must be used for bound- 
ary correspondence on three sides of the map. However, 



for a restricted FOV (such as 90°), we have the freedom 
to change the dimensions of the virtual domain so that 
M = 1 and the QCM reduces to the CM [7]. The material 
parameters of the CM are nearly identical to those shown 
in figure [] by eye, but aberrations in[]are visibly reduced. 
The only remaining aberrations are caused by slight re- 
fraction at the top boundary of the map. Since there is a 
small gradient in the material parameters across the top 
boundary, collimated rays are refracted at different an- 
gles upon entering the transformed domain. Since these 
rays are no longer parallel, they are focused to different 
points on the focal plane. These aberrations are only re- 
duced by extending the boundaries of the map, thereby 
increasing the size of the device. By changing the mod- 
ule in this fashion, we have also increased the material 
discontinuities on the sides of the map. This restricts the 
FOV of the lens. 

We have shown that the QCM technique can be ex- 
tended to the design of rotationally symmetric TO de- 
vices that can be realized using simplified, uniaxial me- 
dia. The Luneburg mapping illustrates several important 
features that are common to the 3D QCM solutions and 
could allow for broadband operation. The dielectric re- 
sponse is mostly positive and moderate everywhere, with 
e'p^^ ranging from 1 to approximately 3.27 and e^ ranging 
from approximately 0.22 to 3.35. However, throughout 
most of the transformed region, e'^ > 1 and we expect 
that we can retain most of the FOV of the lens while ne- 
glecting the portion with values falling below the vacuum 
permittivity. Though magnetic response is required, /i^ 
ranges from approximately 0.15 to one and can thus po- 
tentially be broadband. Broadband diamagnetism has 
been previously suggested, with one possible implemen- 
tation consisting of stacked metal plates oriented to pro- 
vide a magnetic response in the azimuthal directionj20]. 
The tremendous simplification in material properties af- 
forded by the 3D QCM approach increases the viability 
of many potential TO devices. 
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